The Influence of Geometry on the Vibronic Spectra of Quantum Aggregates 
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A study is presented of the localisation of excitonic states on extended molecular aggregates 
composed of identical monomers arising, not from disorder due to statistical energy shifts of the 
monomers, induced by environmental interactions (Anderson localisation), but rather simply due to 
changes in the orientation and geometrical arrangement of the transition dipoles. It is shown further 
that such small changes nevertheless can have a drastic effect on the shape of the vibronic spectrum 
of the aggregate. The vibronic spectra are calculated using the "coherent exciton scattering" (CES) 
approximation whose derivation we generalise to be applicable to aggregates of arbitrary size and 
geometry. 



I. INTRODUCTION 

The quantum behaviour of extended aggregates of 
atomic and molecular monomers, containing from just a 
few, up to thousands of sub-units, is attracting increasing 
attention in chemistry and physics. Prominent examples 
are aggregates of large dye molecules [IHl] , chromophorc 
assemblies describing the photosynthetic unit [3, Q or as- 
semblies of ultra-cold atoms Ja-Q- Clearly, since the in- 
teracting transition dipoles are three-dimensional vectors 
(the interaction is of tensor form) and since the probing 
light can be polarised, the precise geometry of the aggre- 
gate is of crucial importance. These two effects, the rel- 
ative orientation of neighbouring dipoles and their orien- 
tation with respect to the light polarisation, decide both 
the energy dependence of the aggregate density of states 
and the distribution of oscillator strength. These are the 
key elements in deciding the response of the aggregate 
to incident light and how the electronic excitation prop- 
agates subsequently. In the case of aggregates of large 
molecules this has been appreciated for a very long time 
and measured spectra have been used to help infer the 
usually unknown geometry of the aggregate under ob- 
servation. For example, whether dye aggregates exhibit 
a J-band, an H-band, or both, in absorption provides 
strong clues as to aggregate geometry. 

Recently two striking new effects of the influence of 
geometry on spectral absorption have been observed. 
When pressure is applied to cylindrical aggregates of dye 
molecules, one observes significant spectral shifts [l()lT12f 
which are interpreted [ill [l2j as arising from a collapse of 
the cylinder cross-section from circular to elliptical form. 
Similarly in spectral studies of the light-harvesting com- 
plex [13l4l6| it has been shown that departures from a 
circular arrangement lead to strong spectral changes and 
detailed comparison with experiment allows one to pre- 
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diet one particular geometrical arrangement of the chro- 
mophores as most likely. 

Several methods have been put forward to account 
for the continuous broadening of observed spectra due 
to interaction with vibrational and environment degrees 
of freedom. By far the most popular method is to ig- 
nore intra-monomer vibrations and the specific nature 
of external degrees of freedom coupling to the exciton 
and simply to assume that the monomer transition ener- 
gies (diagonal disorder) and/or the inter- monomer elec- 
tronic interaction strengths (off-diagonal disorder) are 
randomly distributed according to some prescription. In 
this model, transition energies with random deviations 
from the mean act effectively as 'impurities' in the ag- 
gregate leading to 'Anderson' localised states which are 
instrumental in inhibiting energy transfer, as seen for ex- 
ample in [13, 03- 

One of the main points to be emphasised below is that 
geometrical conformation, even in a model of monomers 
with identical transition energies, can lead to similar ef- 
fects without random disorder in the couplings. That is, 
abrupt local changes in conformation, e.g. simple bend- 
ing of a polymer chain, give rise to changes in absorp- 
tion spectra and occurrence of localised states. This is 
in contrast to 'Anderson' localisation arising from statis- 
tical fluctuations. We feel that insufficient attention has 
been accorded to the strong spectral effect of such simple 
'geometrical' localisation. Although 'topological' locali- 
sation has been well-studied in other contexts, again it 
has been from a statistical point of view. 

Of course, the formation of localised states, splitting 
off from delocaliscd band of states, due to isolated impu- 
rities is a well-known phenomenon. It occurs not only in 
exciton bands but also in conduction bands and phonon 
bands. However, in our case there is no real impurity, 
just a change in interaction due to changed geometry. As 
illustrative examples of this we will consider the very sim- 
plest cases of a single bend in a linear chain of monomers 
or an elliptical deformation of a circle. 

Many related studies have been presented concentrat- 
ing on statistical distributions of " impurities" . Dynami- 
cal diagonal disorder (i.e. time-dependent fluctuations of 
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monomer transition energies) on circular aggregates has 
been studied by Bakalis et. al. and by Wubs and 
Knoester [20j |. Similarly Freiberg et. al. [2l| studied dis- 
order in the LH1 and LH2 photosynthetic aggregates and 
Barvik et. al [22j calculated lineshapes for rings of BChl 
monomers including statistical disorder in transition en- 
ergies and coupling strengths but only for dichotomic dis- 
order and for averaged nearest-neighbour coupling. In 
a further advance, Fidder et. al. [2!|, the effect of off- 
diagonal disorder in coupling strengths was included via 
a statistical distribution of inter-monomer spacings. Wu 
and Small [HJ also considered both diagonal and off- 
diagonal statistical disorder on circular aggregates with 
reference to the LH2 photosynthetic complex. The ques- 
tion of the influence of elliptical deformation of circular 
structures encountered in LH2 photosynthesis was con- 
sidered by Matsushita et. al. [25| for purely electronic 
transitions and with inclusion of statistical broadening 
by Warns et. al. [Hj]. 

Here wc seek to include vibronic coupling more di- 
rectly than in the statistical models, which essentially 
ignore internal vibrational structure. Even though one is 
dealing with large organic molecules, in statistical mod- 
els they are treated as structureless entities characterised 
only by a given electronic transition energy. As in our 
previous work on spectra and energy transfer [27T - I291 ] . 
we consider first the purely electronic model of iden- 
tical coupled monomers and then the realistic case of 
broad structured spectra arising from interaction with 
both intra-molccular vibrational modes and with exter- 
nal vibrational modes of the environment. We use the 
"coherent exciton scattering" (CES) approximation gen- 
eralised to be applicable to aggregates of arbitrary size 
and geometry. This approach can also take the influ- 
ence of disorder and and temperature implicitly into ac- 
count [30| . In a previous publication (3l| we have tested 
the CES approximation by diagonahsing Hamiltonians 
of the type considered here. For red-shifted (J-band) 
aggregate spectra the CES calculation gives very good 
agreement with exact results for all coupling strengths 
between monomers. For blue shifts (H-band spectra) it 
gives good agreement for weak and strong coupling. Only 
in the intermediate coupling case is the agreement not so 
good, although the main qualitative features of the full 
diagonalised results are still obtained. This theoretical 
study has been performed for the extreme case of a sin- 
gle vibrational mode. Comparison to experiment (where 
the vibrations form a continuum) has shown that in such 
realistic situations the CES approximation is able to re- 
produce the J-band [32j and also the H-band 27 1 in the 
intermediate coupling regime with remarkable accuracy. 

Without the influence of vibrations, the overall spec- 
tral structure is decided by the density and oscillator 
strength distribution of electronic transitions. Hence, we 
will show how the geometry of transition dipole orienta- 
tion and isolated geometrical irregularities can seriously 
influence the purely electronic spectrum of molecular ag- 
gregates. In particular we point out that the orientation 



of transition dipolcs plays a decisive role in deciding the 
characteristics of absorption spectra. Also, depending on 
the orientation, sudden changes in conformation due to 
bends or kinks in regular structures lead to local changes 
in off-diagonal coupling and thus to localisation of wave- 
functions on these sites. In some cases, these localised 
states can still be prominent in the spectrum when vi- 
brational structure is taken into account explicitly. 

The plan of the paper is as follows. In section II we 
derive expressions for the aggregate absorption spectrum 
for the purely electronic case and when coupling to vi- 
brations are included. The simplest case of a bent linear 
chain is considered in section III and that of a circle de- 
formed into an ellipse in section IV before conclusions are 
discussed in section V. 



II. THE AGGREGATE SPECTRUM 
A. the purely electronic case 

The purely electronic model in which the inter- 
monomer interaction is taken to be that of point dipolcs 
located on each of TV monomers is particularly simple. 
The total aggregate Hamiltonian is H = Ho + V, where 
Ho is the sum of N identical monomer electronic Hamil- 
tonians and V is the dipole-dipolc coupling operator. 
First the Hamiltonian is expressed in the localised ba- 
sis of aggregate states, 
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(i) 



where \<f>%) is the excited electronic state and |<^) the 
ground electronic state of monomer n. In this basis the 
Hamiltonian reads, 



H = ^2 e n \ir n )(n n \ + ^ Vnm\ir n ) (n r , 



(2) 



where e n is the transition energy of monomer n and 
the aggregate ground state energy is taken to be zero. 
This is the starting point for statistical models in which 
e n and V nm are considered randomly distributed to in- 
clude vibrational and other broadening phenomcnologi- 
cally. Here, starting first with the purely electronic prob- 
lem for identical monomers, we take all e n to be the same 
and the V nm to be determined from the dipolc-dipole in- 
teraction 
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where jl n is the transition dipole on monomer n and X nm 
is the vector separation of monomers n and to. 
Below the absorption spectrum is expressed in terms of 
the energy-dependent Green function defined by the ag- 
gregate Hamiltonian as, 



G(E) = {E -IL + iS)- 1 , 



0-. 



(4) 
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Similarly, for the non-interacting monomers we have, 

g(E) = (E-H Q +i5)-\ (5) 
The two Green operators are connected by the equation, 
G(E)=g(E)+g(E)VG(E). (6) 
In the localised electronic basis this equation reads 

Gnra — Qn&nm ~t~ Qn ^ ^nn'^n'm* (7) 
n' 

As alternative, which as we will see exhibits geometric ef- 
fects more clearly, we transform to the delocaliscd exciton 
basis, 



\k) 



a k „\ir n ), k = l...N 



(8) 



where the coefficients a kn are determined by diagonal- 
ising the aggregate Hamiltonian H. Note that in this 
paper we use exclusively the subscripts n, m to denote 
the "site" basis and the subscripts fc, k' to denote the ex- 
citon basis. In the exciton basis one readily sees that the 
Hamiltonian has matrix element 



|H|fc') = (e + C k )S. 



kk 1 



where, 



k a mk 



(9) 



(10) 



In the purely electronic, one-level model one has the 
simple form g n (E) = (E — e n + i(5 + ) _1 . Hence, for identi- 
cal monomers, Eq.© simplifies since g n — g is indepen- 



dent of n, i.e., 



(11) 



When transformed to the exciton basis this equation 
reads. 



Gkk' = gSkk' + gCkGkk' 



Gkk' — 



a-gCk) 



Skk> 



(12) 



(13) 



Since g(E) = (E — e + i5 + ) 1 , from the above equation 
(TTBl one has 



Gkk'(E) 



1 



(E - e - C k + iS+) 



6 kk >, (14) 



as must be since H, and therefore G, is diagonal in the 
exciton basis with eigenergics e + Ck ■ Although trivial in 



this case, the above formula ((13)) will be useful when vi- 
brations are included. The shape of the absorption spec- 
trum is given b y — ImA(E) where the spectral function 

A(E) is hi mil, 



A(E) = J2(^n)(P m ^)G nm (E) 



(15) 



where fl n is the transition dipole on monomer nand e is 
the light polarization vector. Transformed to the diago- 
nal exciton representation, this equation becomes, 



A(E)=J2\(e-tk)\ 2 G kk (E)=J2\ 



9(E) 



(i-g(E)CkY 

(16) 

where the effective transition dipole in exciton state k 
is given by fl k = X) ti M»i a «.r This equation makes clear 
how the geometry influences the absorption spectrum via 
both fl k and in the interaction C k . Clearly, if the intra- 
monomer interaction is switched off, then the monomer 
spectral function is 



M(E) = (J2\(e-M\ 2 )g(E) 

n 

or, equivalently, in the exciton basis, 

M(E) = (Y,\(e-tk)\ 2 )g(E) 



(17) 



(18) 



B. Inclusion of Vibronic Coupling 

To calculate realistic spectra one must include the 
broadening due to coupling to both IM and EM vibra- 
tional degrees of freedom. This requires inclusion of these 
degrees of freedom in the definition of the non-dipole cou- 
pled monomer Hamiltonian Ho. When this is done, the 
Green operators g(E) and G(E) become dependent on 
these degrees of freedom. Then the spectral function (fT"5)) 
must be generalised to [H, 0, HU 



A(E) = J2(e-M(fI m ^)(G nm (E)), 



(19) 



where (■■•■) denotes an integration over the full vibra- 
tional ground state of the aggregate. Note however, that 
Eq. (fl6]) cannot be so simply generalised since G is no 
longer diagonal in the pure exciton basis. Nevertheless, 
we can still transform to this basis, when Eq.(|16j) is re- 
placed by, 



A{E) = Y,{Z-tfkWk.?)(G k k-{E)) 



kk' 



Similarly in this basis, Eq.([B]) becomes 

(G kk ') = (9kk>) + ((sVG) kk ,) 



(20) 



(21) 
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FIG. 1: The energy eigenvalues (left column) and electronic spectra (middle column) for a linear chain of monomers with 
/t = (0,0,1) and bend angle $ = 0°, 120°, 135° in a), b), c) respectively. The vibronic CES spectra (right column) are 
calculated with V = 150 cm" _1 (red dotted), V = 300 cm -1 ( green dash-dotted) and V = 450 cm 1 (blue dashed). 



In general this equation cannot be simplified further. To 
proceed, two strategies are possible. One is to work fully 
numerically and simply expand the Hilbert space to a 
basis which is not purely electronic but which includes 
all vibronic degrees of freedom. Then the Hamiltonian 
and G(E) again are diagonalised in this expanded ba- 
sis. In view of the explosion in the number of vibronic 
basis states however, this strategy is only possible for 
very small aggregates and only a few vibrational modes 
[3ll. I3614421 ]. Alternatively one can seek approximations. 
In previous work we have shown that the 'coherent exci- 
ton scattering' (CES) approximation gives agood repro- 
duction of experimental aggregate spectra |27l. r28l |32|. l3o| 
and agrees with exact diagonalisation results for a variety 
of examples (30l . l3l| . In this approximation, the operator 
g in the second term on the right hand side of Eq. (f2"Tj) 
is replaced by its ground-state vibrational average (g). 
Then, if again we assume that (g n ) is independent of n 
due to identity of the monomers, it is easy to show that 
Eq. (|2"Tj) reduces to, 



(G kk >) = (g)5 kk , + (g)C k (G kk ,) 



(22) 



with Cfc given by Eq. (|TU)) . Then one has again the simple 
result, 



(GW) 



(9) 



(1 - (g)c k ) 



and 



A(E) = ^2\(e.p k )\ 2 - 



(9(E)) 



(l-(g(E))C k y 



(23) 



(24) 



which arc identical to the purely electronic case Eq. (fT3")) 
and Eq. (|TB|) except that g is replaced by its vi- 



bronic counterpart (g). Similarly, the spectrum of non- 
interacting monomers is given by Eq. (|17|) with g re- 
placed by (g). This has the important consequence that 
all monomer and aggregate spectra are broadened due 
to simultaneous vibrational excitation and environment 
coupling. The input to the CES spectral calculation 
is the complex function (g(E)) which we will take to 
be of such a form as to reproduce typical experimen- 
tal monomer spectra composed of several, continuously- 
broadened peaks as shown in Fig. ([T]) right column, black 
curve. 



III. THE BENT LINEAR CHAIN 

First of all we consider a rigid linear chain of equally- 
spaced identical monomers. A bend is introduced at one 
monomer (we take number 12) to give a V-shaped chain 
in which the two half-chains are at a mutual angle $ 
(see the insets on Fig.([T]), where we have taken TV = 19 
monomers and located the bend on monomer 12). Se- 
quences of such bent chain segments are a good model 
for certain types of conjugated polymers, as studied for 
example in Ref. [43[. Although we consider interactions 
between all monomers on the chain, it is clear that the 
interactions between the monomer at the vertex and its 
nearest neighbours will suffer the largest change upon 
bending. Indeed we will find that this simple symme- 
try breaking causes the appearance of localised states 
centred around the vertex. The linear case $ = is 
shown in Fig. (fjji). Here the dipoles have been taken 
to be aligned perpendicular to the chain. The nearest- 
neighbour interaction from Eq. Q for this geometry will 
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be taken as the unit of energy in all cases. Restricting to 
nearest-neighbours alone, in this unit the exciton band 
would occupy the domain —2 to +2 . Including all cou- 
plings shifts the band slightly, as seen in Fig. (QJ,) and 
discussed e.g. in j23j . 

For each geometry we calculate the excitonic eigenval- 
ues and eigenvectors (wave-functions) by direct numeri- 
cal diagonalisation and in the first and second columns of 
Fig. (JTJ) we present the energy eigenvalues and oscillator 
strength 'stick' spectra respectively. Here and in subse- 
quent figures showing spectra, the undistorted nearest- 
neighbour value of V n , n +i hi Eq. (fTU)) , which we call the 
coupling strength V, is taken as the energy unit for the 
stick spectra in the second column. Note also that we 




12 17 2 12 17 



Monomer number m Monomer number m 



FIG. 2: Wavefunctions, left column, and modulus squared, 
right column, of a bent linear chain with (f> = 135° as a func- 
tion of monomer number. The energy of the states is increas- 
ing from bottom to top. 



have normalised the oscillator strengths to sum to unity 
in each case. Only in assembling the CES vibronic spec- 
trum shown in the third column do we need absolute 
values. Here we have taken as typical monomer spec- 
trum that of pinacyanol, stretching from roughly 17000 
to 20000 cm -1 and calculated the aggregate spectra using 
Eq. (j2"4")l . Then the couplings Ck must also be absolute, 
which is achieved by fixing V. The realistic size of this is 
determined by the width A rj 2000 cm -1 of the monomer 
spectrum. In Fig. ((T|) and in subsequent figures we plot 
the continuous vibronic spectra in absolute energy units 
of cm -1 . 

For the straight chain, shown in the top row of Fig. [TJi, 
the results are well-known. The energy spectrum is a 
smooth function of k and one level at the band edge (a 
blue-shifted so-called H band in this geometry) carries 
almost all the oscillator strength. Inspection of the wave- 
functions (not shown) shows that all levels are delocalised 
over the whole chain. In the right column we show the 
CES spectrum for the monomer and for weak, (V = 150 
cm -1 ), intermediate (V — 300 cm" 1 ) and strong (V = 
450 cm -1 ) coupling. As the coupling strength increases, 
the monomer spectrum metamorphoses into a typical, 
relatively broad blue-shifted H band. Note, that the dis- 
tinction weak / intermediate/strong is not based on the 
Simpson-Peterson criterion j44j but is just introduced to 
denote the relative strengths of the couplings. 

We take the bent chain to lie in the (x, y) plane. Cal- 
culations for finite $ show little change both in the pure 
electronic and the vibronic spectrum so long as $ < 90°. 
However, for obtuse angles there begin changes in the 
spectrum, as shown in Fig.(Q}3) for $ = 120°. Here the 
highest level carrying most oscillator strength has begun 
to separate from the rest of the exciton band levels. This 
is seen very clearly in Fig. ([TJ:) for $ = 135°. Now two 
levels, one above and one below the band, have split off 
. An inspection of the wave-functions (Fig. ^ allows a 
simple explanation in that the two levels which split off 
are strongly localised on only two monomers. This is 
our first example of 'geometrical' localisation. The sharp 
bend in the chain gives a strong interaction between the 
two monomers adjacent to the vertex monomer. This 
leads effectively to the formation of an isolated dimer 'im- 
purity' between these two monomers. The pair of states 
split off from the band in Fig. (fib) are just the symmet- 
ric '+' and antisymmetric '— ' states of this dimer whose 
wave-functions are localised on just these two monomers, 
see the lowest and highest energy states of FigJ2] For this 
dipole orientation, only the high-energy symmetric state 
carries oscillator strength, as seen in Fig.([TJ;). The ac- 
tual vertex monomer has now become effectively isolated 
from the rest of the chain and leads to a further state 
localised on just this monomer. However this state sits 
in the middle of the band (state number 8) and carries 
little oscillator strength. All other states are delocalised 
but, since the chain now consists of two halves, they are 
predominantly confined to one or other half of the chain, 
as shown clearly in Fig. (2). Apart from the highest en- 
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a) 




16000 18000 20000 22000 



b) 




state numbc 




energy 



16000 18000 20000 22000 

energy in cm - 1 



FIG. 3: Eigenvalues and spectra as in Fig.(l) for a bent chain with </> = 135° and dipoles with orientation a) jl — (0, 1, l)/y/2 
and b) fi — (1, 0, 0). In the right column the CES spectra are for the same values V of the coupling as in Fig.(l). 



crgy dimer '+' state the delocalised in-phase second and 
third states from the top of the band are the only ones 
to carry oscillator strength, as seen in the middle col- 
umn of Fig. (jTJj) . Inclusion of vibrations, right column of 
Fig. (JTfc) shows the expected formation of an H band for 
weak and intermediate coupling, but for strong coupling 
the isolated, dimer peak gives rise to a more pronounced 
shoulder on the spectrum. Note however, that due to the 
inclusion of vibrations the highest electronic absorption 
peak is now accompanied by a complicated vibrational 
structure. This in turn makes it difficult to distinguish 
electronic and vibrational contibutions from measured 
spectra, similar to the situation in dimers (45l . 146^ . 

To illustrate the spectral changes resulting from rela- 
tively simple changes in dipole orientation we now con- 
sider a bend in which the dipoles are still perpendicular 
to the chain axis but now lie in the plane of the V-shaped 
bent chain, not perpendicular to it as above. The first 
case we consider is where the dipoles have [i perpendic- 
ular to the chain but with equal components in- and out 
of-planc, shown in Fig. [3^, for $ = 135°. Of course the 
straight chain spectrum is again that of Fig. and as 
the chain is bent, up to $ « 120°, the situation is roughly 
as in Fig.lHJ)) with the formation of a broad H band for 
strong coupling. These absorbing states are delocalised. 
However, for <E> = 135° the split-off dimer localised states 
again become evident above and below the band, but 
now both symmetric and antisymmetric localised states 
carry oscillator strength. Then as shown in Fig. (JTJ, for 
<I> = 135° sharp peaks appear both above and below the 
exciton band. In fact the sharp peak below the band 
(appearing for strong coupling) is reminiscent of a col- 
lective fully delocalised J-band exciton state but we em- 
phasise that it is of completely different character in that 
its wave-function is localised on just the two monomers 
adjacent to the vertex. 

This similarity of spectral appearance of localised and 
delocalised states is also seen when we consider a linear 
aggregate showing a J-band and not an H-band as in 



Figs. fl]). This is the case when the monomer transition 
dipoles are all aligned initially parallel to the straight 
chain as shown in Fig. ^p). The electronic spectrum is 
dominated by the lowest energy level where all dipoles are 
in phase. In the vibronic spectrum of the straight chain 
(not shown) a clear J-band based on this state forms as 
the coupling strength increases. Again, for bend angles 
up to « 120° we find similar spectra to zero degrees i.e. 
this single peak dominates. However, at larger angles, 
(the case <I> = 135° is shown in Fig. (J3Jd) as example), 
again localized dimer states appear above and below the 
band but now the absorbing '+' state is located below the 
band. Also, since the strong vertex interaction effectively 
splits the chain into two linear segments, examination of 
the wave-function shows that the J band is due to two 
collective states (the second and third states in the stick 
spectrum of Fig. (JSJd) again each largely confined to a 
separate half of the chain. This gives a J-band in strong 
coupling. On the low energy side of these two peaks ap- 
pears the localised symmetric dimer state. Note that in 
the vibronic spectrum of Fig. (|3p) these features are only 
evident when the coupling is strong enough to move their 
energies below the region of monomer absorption, when 
the localized state appears as a low-energy shoulder on 
the J-band. 

The message of the preceding section is that even in the 
very simple case of a single fold in a linear chain new 
spectral features appear which depend crucially on dipole 
orientation and which make interpretation of spectral fea- 
tures in terms of localised and delocalised states difficult. 
This will be more so in the case of the complicated geome- 
tries occurring in practical cases in biological aggregates 
for example. 
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IV. CIRCULAR AND ELLIPTICAL 
GEOMETRIES 

Here we will emphasise the effects of dipolc orienta- 
tion. The first aim is to show that even in completely 
symmetric geometries, a changing dipolc orientation can 
affect the monomer interactions dramatically and lead 
to corresponding dramatic changes in the energy loca- 
tion of eigenstates and their light absorbing and emitting 
properties. This also will show that in general it is abso- 
lutely necessary to take all interactions into account: in 
some cases the much-used nearest-neighbour approxima- 
tion fails completely. 

The simplest non-linear geometry (but which is practi- 
cally important in photosynthesis) is a circular arrange- 
ment of dipoles. As additional simplification, first we 
consider equal spacing of dipoles all lying in the plane 
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FIG. 4: Eigenvalues and stick spectra for a circular aggregate 
with dipoles lying in the plane at angle <f> w.r.t. the tangent, 
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of the circle. Then we examine the energy and absorp- 
tion spectrum as a function of the angle <f> of the dipolc 
w.r.t. the tangent to the circle i.e. = 0° is tangen- 
tial and (j> = 90° is radial. The results are shown in 
Fig. flU). Each pair of figures is for a different 4> orienta- 
tion and shows in the left column the geometry and the 
exciton energy) level distribution (w.r.t. e = 0) between 
the lower and upper band edges. In the right column is 
shown the corresponding absorption stick spectrum. The 
results can be understood analytically from the solution 
using nearest-neighbour coupling only. In this geometry 
the exciton states are all delocalised over the whole ring 
and the eigenenergies are (for ./V even) 



N/2 



E k = 2 ^2 cos(fcn)K,i + (-1)^/2+1,1, 



with 



n=2 



I 277 ■ 
K = ?, 

N 



3 =0,1 A 1. 



(25) 



(26) 



As seen from Fig. (J4^i) for = 0°, the levels are doubly 
degenerate apart from the extreme j = (all monomers 
in phase) and j = N/2 (successive monomers out of 
phase) levels which lie at the bottom and top of the band 
respectively. Only the two degenerate levels j ' = 1 and 
j = N — 1 carry oscillator strength. In particular the 
lowest j = level in Fig.(|2Ji) is forbidden. Clearly the 
presence of this dark state below the strongly-absorbing 
J-band states is important for the fluorescence properties 
of aggregates with this geometry. For <fi = 90° the energy 
level ordering is simply reversed and the absorbing states 
form an H-band with the dark j = state lying above 
them in energy. 

As (j) increases from zero to approach the 'magic' an- 
gle of ~ 54° around which nearest-neighbour coupling 
vanishes, the major spectral change is that the lowest 
non-degenerate j ' = dark level moves up in energy. In 
Fig.(|3]3, first column) for <fr — 48° it is state number 3 
and for Fig.flU:) <fr = 53° one sees this non-degenerate 
level as state number 9. Note that in this region the 
nearest-neighbour interaction goes through zero and the 
higher-order interactions take over. By <j> = 54° the j = 
level has reached the top of the band (Fig.(j4ji)) and is 
being followed by the allowed levels which are now lo- 
cated in the centre of the extremely narrow band. This 
band reversal is extraordinarily sensitive to orientation 
and already by (j) = 55° (Fig.(|4^)) the H-band config- 
uration has been reached and there is little qualitative 
change up to = 90° of Fig.(|3F). We do not show the 
broadened CES spectra but it is clear that they change 
from sharp, fully delocalised J-band character for tan- 
gential orientation, through almost unchanged monomer 
character around 54° orientation to H-band character for 
radially oriented dipoles. 

Having illustrated the extreme sensitivity of spectra 
to dipolc orientation on a perfectly symmetric circular 
structure, we examine next the effect of symmetry- 
breaking by deforming the circle successively into 



8 




FIG. 5: As in Fig. [4]the dipoles are in the ring plane and (f> — 0° wrt the tangent on the ring. The upper row shows / = (i.e. 
a ring) and the lower row a strong deformation / = 0.7. The CES spectra in the right column are for V = 50, 100 and 150 
mf 1 (red, green , blue). 



an ellipse of increasing eccentricity. Rather than ec- 
centricity, we will use the simpler 'flattening factor' 
/ = 1 — (b/a), where b and a are the minor and major 
axes of the ellipse, respectively. Hence, / = for a circle 
and / = 1 when the circle is squashed flat. We consider 
two dipolc arrangements on the circle, the first of which 
is relatively insensitive to deformation, the second of 
which shows strong sensitivity. 

The first geometry is the tangential arrangement of 
Fig(j4^,), shown again in Fig.([5K). In the right-hand 
column we now show the CES vibronic spectra. For 
a circle, the spectrum shows the characteristic J-band 
structure for strong coupling. In fact as the circle is 
squashed into an ellipse there is almost no perceptible 
change in the CES spectrum. Therefore, in Fig. we 
show only the rather extreme case / = 0.7, where the 
similarity to the circular case is confirmed. The reason 
for this can be explained by examination of the exciton 
wave- functions (not shown). As the ellipse is deformed, 
states localised at the corners of the ellipse are formed 
but they remain in the middle of the exciton band and 
do not absorb appreciably. The two long segments of the 
ellipse become more and more like two linear chains. In 
this orientation, these have the same J-band absorption 
characteristics as the circle. The small changes that do 
occur are washed out in the CES spectrum, so that the 
right hand columns of Figs.([5£i and \5]p) show spectra 
which arc quite similar. 

As an arrangement which is extremely sensitive to 
symmetry-breaking we choose the case of in-plane 
dipoles with 4> = 54° orientation of Fig.Qd. This is 
reproduced as Fig.^). Here the nearest- neighbour 
coupling is essentially zero. One notes the two allowed 
levels sitting in the middle of a narrow band with the 
non-absorbing ,7 = level split off at the top of the 
band. The CES aggregate spectrum (right column) is 



then almost the same as the monomer spectrum. 
As the circle is flattened, spectacular changes occur as 
the nearest-neighbour coupling switches on. Already for 
/ = 0.2, Fig.([6p), the two degenerate allowed levels have 
split to be the bottom and, apart from j = 0, top levels 
of the band. These states are still fully delocalised on 
the ellipse. However, as flattening progresses Figs.([Bb,c), 
three new features arise. Firstly, the allowed levels 
acquire a degenerate partner and become the lowest 
and highest energy levels. Secondly, they become 
localised at the corners of the ellipse. Inspection of the 
wave-functions shows that, as in the linear case, they 
are nothing other than the '+' and '— ' dimer levels, now 
doubly-degenerate and both dipole-allowed because of 
the orientation. The rest of the ellipse behaves as two 
linear segments so that the third feature to appear is a 
central exciton band with a new, fully delocalised J-band 
state at its lower edge and a weak H-band absorption at 
the upper edge. These features of dimer-pair formation 
and the splitting off of localised states to leave a central 
band, are seen clearly in the first and second columns 
of Fig. (jBji). The CES spectrum in the third column 
shows, in strong coupling, an apparent double J-band 
structure and a broad blue-shifted H band. However, 
we emphasise the completely different characters of the 
two apparently similar J bands. The lower is a localised 
state due to geometric disorder, the other is a delocalised 
exciton state. Comparison of the CES spectra of Fig. ((St) 
and Fig.JSb), both for / = 0.7, emphasises the strong 
changes that can occur simply due to dipole orientation. 

As a last example we consider out-of-plane orientation 
of the dipoles, with a z-axis perpendicular to the ellipti- 
cal plane. We take all three (x, y, z) dipole components 
to be equal corresponding to <j) — 45° and polar angle 
= 55°. From Fig. ([7^,) one sees that the circular ge- 
ometry energy-level distribution (first column) is similar 
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FIG. 6: As in Fig. ((5} but for in-plane dipoles with orientation (j> = 54° w.r.t. the tangent and aggregates with a) / = 0, b) 
/ = 0.2, c) / = 0.4, d) / = 0.7. The CES spectra in the right column are for V = 150, 300 and 450 cm -1 (red, green , blue). 



to the <j) = 54° case of Fig. (JS^i), but the stick spectrum 
of the second column is quite different. Now the two 
j = 1 absorbing levels with polarisation in the plane ap- 
pear at the high-energy H-band side. In addition, the 
fully-symmetric, fully-delocalised J = level now carries 
oscillator strength polarised perpendicular to the plane. 



Upon flattening to an ellipse, the changes mirror those 
of Fig.©, except that the j = level maintains its 
dipole-allowed absorption throughout, since there are no 
changes in geometry in the perpendicular direction. For 
/ = 0.4 the j = 1 levels have acquired degenerate part- 
ners and separated, as states localised at the ellipse edges, 
from the bottom and top of the band of remaining levels. 
Of these levels, the lowest is again a delocaliscd exci- 
ton state carrying most of the oscillator strength and the 
highest levels now carry oscillator strength mostly po- 
larised in the perpendicular direction. The result for the 
strong-coupling vibronic spectrum Fig. ([7]:, right column) 
shows again a J-band structure but now with a very pro- 
nounced blue-shifted H band apparent. 



V. RELEVANCE FOR EXPERIMENTS 

This study was motivated by observations of strong 
spectral changes accompanying pressure change giving 
rise to a presumed flattening of cylindrical aggregates 
composed of different substitucnts of 5,5', 6J3' - tetra- 
chlorobenzimidacarbocyanine (TDBC) dyes [lllQjl- The 
following characteristics under increasing pressure were 
noted; 

a) the non-aggregated monomer Stokes' shift increases 
linearly with pressure. 

b) Unlike other studies on cyanine dyes (ic| and por- 
phyrins [47[ , the red shift of the J-band of the aggregate 
is not monotonically linear due to compression of the 
inter-monomer distance. Rather a linear shift is followed 
by saturation as pressure increases. 

c) Nevertheless the aggregate Stokes' shift does in- 
crease monotonically with pressure. 

d) As pressure increases a new absorption band arises 
as a shoulder on the spectrum at the low-energy side of 
the J-band. 

Spitz and Daehne [l2[ suggested that these changes 
arise from a collapse of the cylinder at a given pressure 
and formation of a localised state at the 'kink' in the 
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FIG. 7: As in Fig. (6) but for out-of-plane dipoles having <f> = 45° and polar angle 6 = 55° in the circular aggregate. The 
figures are for a) / = 0, b) / = 0.2, c) / = 0.4. The CES spectra are with V = 750 cm" 1 , 1500 cm -1 and 2250 cm" 1 (red, 
green , blue). 



now almost elliptical cross-section. The localisation of 
excitation on just two monomers implies that the Stokes' 
shift behaves in the same way as for the monomers 
(J-band emission shows a small Stokes' shift) and also 
can account for the new band emerging below the J- 
band. Unfortunately, due to the lack of knowledge of the 
precise geometry and particularly dipole orientation in 
such aggregates, we do not consider it realistic to present 
our detailed calculations showing these characteristics. 
Suffice it to say that the calculated spectra of elliptical 
conformations such as shown in Fig. (H5J: and d) where 
a dipole-allowed localised state arises below the J-band 
level do lend qualitative support to the interpretation of 
Spitz and Daehne. 



VI. CONCLUSIONS 

We have examined the influence of dipolar orientation 
and symmetry-breaking on the energy and light absorp- 
tion spectra of very simple linear and circular aggregates 
of identical monomers. The major features to appear 



are a strong dependence of absorption characteristics 
on orientation and the isolation of states localised 
around bends in regular structures. In many cases these 
localised states are split off from the top and bottom 
of the exciton band and may carry significant oscillator 
strength. Clearly the presence of such states, either 
dark or absorbing, below the exciton band edge can 
seriously affect the fluorescence properties and Stokes 
shift of the aggregate. The flattening of a circle also 
leads to the formation of localised states around the 
ellipse extremities. Depending upon dipole orientation 
these can lead to considerable broadening and blue shift 
of the spectrum. The results shown are admittedly 
model calculations, although we feel that the inclusion 
of vibrations in the CES approximation is realistic since 
in previous work we have obtained very good agreement 
with experiment. The aim of this study has been to 
draw attention to the localisation occurring due to very 
simple geometrical deformation of symmetric aggregates 
and to point out that this localisation can lead to new 
spectral structures. 
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